Time-reversal symmetry in nonlinear optics 
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The applicability of time-reversal symmetry to nonlinear optics is discussed, both from macro- 
scopic (Maxwell equations) and microscopic (quantum theoretical) point of view. We find that only 
spatial operations can be applied for the symmetry classification of nonlinear optical processes in 
magnetic, in particular antiferromagnetic, materials. An example is given where both operations 
(time reversal and a spatial operation) can yield different results. 
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Symmetries determine several important properties of 
a crystal, in particular its optical response. In magnetic 
materials, time-reversal is believed to be of fundamen- 
tal importance since this operation reverses all magnetic 
moments [Q-^. However, the consequences of applying 
time-reversal are more profound than a simple inversion 
of localized magnetic moments. As it will turn out in this 
paper there is a deep interrelation between the absence 
of conventional dissipation in even-order (e.g. second) 
harmonic generation and the influence of time-reversal 
on spin ordering. This brings about a subtle difference 
between time-reversal and spatial symmetries in nonlin- 
ear optics. The benefit of this difference makes opti- 
cal second harmonic generation (SHG) a rather unique 
probe of antiferromagnetism, while linear optics (where 
dissipation in the conventional sense is possible) is blind 
for such balanced spin structures. The recent discus- 
sion about the influence of micro-irreversibility on macro- 
reversibility and reciprocity ( |^-^) shows that the is- 
sue of time-reversal, although extensively discussed, is 
far from being understood. 

The theory of nonlinear optics has been developed 
since the 60s. The pioneering work of Armstrong et al. |^] 
describes the propagation of a light wave through a non- 
linear medium, where the energy may be converted from 
the fundamental frequency to higher harmonics (or vice 
versa). An exhaustive description of nonlinear optical 
phenomena is contained in the fundamental books by 
Bloembergen and Shen In these works ( a 
unique flow of time is tacitly assumed, while magnetism 
is entirely absent. Consequently the issue of time re- 
versal is not essential for these authors. The discussion 
of magnetism has been brought to nonlinear optics by 
Pan et al. jl^] and Hiibner et al. In these papers, 

time-reversal was applied to reverse the localized mag- 
netic moments, since the discussion was focused on fer- 
romagnetism. However, the experimental observation of 
antiferromagnetic (AF) domains in Cr203 by Fiebig et al. 
| p2[ and the subsequent theoretical analyses by Muthuku- 
mar et al. |jl^ and by Dahn et al. challenged the va- 
lidity of time-reversal for the symmetry analysis of optical 



processes. Since the inclusion or absence of time-reversal 
in the theoretical analysis of SHG from antiferromagnets 
yields different predictions of the experimental results, 
the issue is shifted from academic interest to practical 
relevance. The importance of the theoretical analysis 
of SHG from antiferromagnets is tremendously growing 
due to the unique capabilities of this method in prob- 
ing buried AF layers, which in turn is important for the 
characterization of recently upcoming magnetoelectronic 
devices such as tunneling magnetoresistive junctions. 

In considering the time-reversibility of an experimental 
situation, three approaches are possible: (i) time-reversal 
is applied to the sample, but all the processes resulting 
from the experiment are unchanged. In particular, the 
magnetic moments in the sample are reversed, but the 
direction of the light propagation through the sample is 
not affected. This approach is presented e.g. in Jl5| , p^ . 
We consider this approach as incomplete, since it does 
not equally treat the sample and the light propagating 
through it. (ii) The second approach, usually encoun- 
tered in the so-called Sagnac-interferometry, addresses 
time-reversal by reversing the propagation of the light 
through the sample (see, e.g. p7|-p^). Clearly, such pro- 
cedure probes the reciprocity of the sample rather than 
its time-reversal symmetry. It can also be proven that 
the second approach is equivalent to the first one. (iii) 
According to the third approach, presented e.g. in po| , 
time reversal acts on both: the sample and the experi- 
mental setup. In this letter, we will follow approach (iii). 

In the processes of even-order harmonic generation, 
dissipation in the conventional sense, converting radia- 
tion into heat, does not exist, since the energy loss of the 
electromagnetic field is the time average 



dt 



(1) 



which vanishes for SHG (and all even-order harmonics), 
since 



P{t) - Poe'"* 



(2) 
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Here, P and E denote the polarization of the medium 
and the electric field, respectively [^. The lack of dis- 
sipation in the conventional sense does not mean that 
the process of SHG is reversible. Already the analysis 
by Armstrong et al. assumes a unique time direction. 
There, the nonlinear polarization P^-^ and the electric 
field E3 of a light beam resulting from Sum Frequency 
Generation at a point tq is given by: 



^j(Ak-ro+A0) i(k3ro-W3<+<^3) 



Re 



3i(k3-ro-W3t+03) 



(3) 



(4) 



see eqs. (3.1) and (3.2) of Ref. 0. Here, ^3 and de- 
scribe the frequency and wave vector of the generated 
light {uj^ = LOi + 102 and ks « ki + k2). The authors 
introduce the idea of "work done on this wave" by the 
nonlinear polarization of the medium, equal to 



W3 
27r 



E3 T^^dt 



cycle 



dt 



= -c^3E3P^^^(cc;3,out-of-phase), (5) 



if the polarization is exactly 90° out of phase with the 
electric field (which requires that Ak^z + A(j> = tt/2). 
The work done on the generated wave determines the 
direction of time. This presents a new kind of dissipa- 
tion, namely "dissipation in the frequency space" , which 
invalidates time-reversal symmetry. 

This fact becomes even more obvious if one takes the 
global picture of SHG. Radiation acting on an ensemble 
of atoms may excite and deexcite them in many ways 
simultaneously. Thus contributions of many frequencies 
are always present (see Fig. |l|(a)). One has a unique 
source of ui light but several detectors for beams of dif- 
ferent frequencies: 2lu, Sw, etc, resulting from sum fre- 
quency generation (in particular SHG); linearly propa- 
gating u> light; and a DC current resulting from difference 
frequency generation. This is due to the expansion of the 
source term (polarization P) in terms of the electric field: 



Pl+P2 + ...= 



(6) 



Imposing time reversal, the detectors become sources and 
vice versa. Thus, in the time reversed process, one ends 
up with a single detector, the one which receives the light 
of frequency uj (Fig. |l|(b)). In order to obtain this single 
frequency one has to redirect all these (previously gener- 
ated) beams back to the sample, conserving their phases. 
The source term now becomes: 



(w): E^'^^E^'^) + x^^\2uj): E^^'^^E 



(2c^) 



(7) 



Since the phases of the now incident electric fields are 
the same as for the previously outgoing electric fields, all 
the terms but those with x^'^'' cancel (which means that 
in the outgoing light one now has only the contribution 
at the frequency lo) and the original situation at the in- 
put of the process is restored. This description, though 
mathematically correct, is physically invalid, since there 
is no practical way to detect an infinite array of frequen- 
cies along with the phases and to revert it with arbitrary 
accuracy (Fig. gc)). Tracing out the "bath" de grees of 
freedom (frequencies other than u and 2a;) causes a tran- 
sition from a pure to a mixed state of the system, which 
means that some memory is lost. This happens because 
the traced subsystem and the bath are not statistically 
independent ||2^. Thus, in any practical situation, there 
is no possibility to generate only the frequency lu out of 
a whole array of frequencies. The process of SHG looks 
different in (-t) than in (t). Such a process is called dy- 
namical. 

As stated before, there is no dissipation in the process 
of SHG in the usual meaning, i.e. the amount of energy 
in the radiative form is constant. However, there is a 
transfer of energy between the frequencies, in particular 
energy flows from the frequency u to other frequencies 
(see Fig. H). We call this dissipation in frequency space, 
in contrast to the more usual dissipation in real time. 
Dissipation in frequency space can mix real and imag- 
inary parts of the nonlinear susceptibility tensor. The 
distinction between these two types of dissipation is of- 
ten encountered in the literature. We consider them here 
on an equal footing stating that the presence of any of 
them (in our case it is the dissipation in frequency space) 
causes the system to have dynamical and thus irreversible 
properties. In this case, time-reversal does not apply to 
the symmetry analysis |ll3| , |l^ , ^ . 

So far we have reasoned that the time-reversal oper- 
ation has to be excluded from the symmetry analysis 
of SHG. However, magnetism may bring an additional 
complication, since the magnetic spin structure is an ad- 
ditional aspect the symmetry analysis must account for, 
and it is the time-reversal which is conveniently applied 
to flip the local magnetic moments. This is, however, 
not correct: it is the classical covering symmetry ]2^ ] 
of the magnetic crystal which should be addressed in a 
symmetry analysis rather than the quantum-mechanical 
symmetry of the wavefunctions p^ . This means that 
the operation applied to reverse the localized magnetic 
moments should be performed in real space rather than 
Hilbert spin space. Consequently, time-reversal cannot 
be used for the symmetry classification of magnetic mo- 
ments. 

Taking into account that time-reversal is not suitable 
for the description of dynamical phenomena, one needs 
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an operation which merely flips the locahzed magnetic 
moments without inverting the time-flow. This can be 
accomphshed by purely spatial point-group operations. 
In many antiferromagnetic crystals a simple translation 
by a lattice vector reverses the magnetic moments. In 
many ferromagnetic and antiferromagnetic systems this 
may be accomplished by a mirror operation. The spa- 
tial operation, which reverses the localized magnetic mo- 
ments, is called by us "moment-reversal" . This operation 
is obviously unitary, in contrast to the time-reversal op- 
eration. Consequently, one does not need to invoke the 
time-reversal operation to describe the full symmetry of 
magnetic crystals. 

Next, we support our reasoning by an example where 
the application of time-reversal and "moment-reversal" 
in the symmetry analysis yields different results (see Fig. 

Let us assume a spin structure with two domains, 
A and B, related to each other by spin- reversal A 
symmetry analysis, similar to the one in pTj , provides 
us with the set of nonvanishing elements of the nonlinear 
susceptibility tensor (i.e. x*-^^ tensor) along with the par- 
ities of these elements. Let us assume that for a certain 
experimental geometry only two tensor elements, called 

(2) (2) 

Xo and Xe , contribute to the resulting SHG hght, and 
that Xo^' is odd while xi^^ is even in the domain opera- 
tion. The intensity of SHG light at a flxed polarization 
is given by: 

Ip-\{xi'^f + ixi'^r±2x^^^-X^^^\ (8) 

where "-I-" stands for domain A, "-" for domain B. In the 
conventional approach, where time-reversal is the oper- 
ation mapping domains into each other, Xo must be 

(2) r~i 

purely imaginary and Xe purely real (Fig. |3|(a)), since 
the operation of time-reversal is anti- unitary pq | . In this 
traditional approach, the flrst two components of the sum 
in eq. (^) are real, while the last one is imaginary. Be- 
cause it is the modulus of the whole sum that determines 
the output intensity, the domain contrast is lost since 

|a-fib| = |a-ib|. (9) 

This is not the case if one uses the spatial operation of 
"moment-reversal" for the symmetry classification, since 
then both tensor elements x'o^ ^iid Xe^'' a-re just complex 
numbers without any constraints on their relative phase, 
see Fig. and domain imaging is possible, as described 
in . Consequently, the symmetry analysis yields very 
different predictions if one uses time- or spin-reversal. In 
the limit far from resonances, however, the phase dif- 

(2) (2) 

ference between Xo and Xe approaches 90°, and the 
domain contrast is lost also in the "moment-reversal" de- 
scription (in agreement with experiment Jl2|). 

Finally we would like to remark on the validity of 
previous work on the group-theoretical classification of 
(magneto-) optical tensors. According to Pan et. al. 



|1C| , the time-reversal operation, because of its anti- 
unitarity, forces the tensor elements to decouple into mu- 
tually exclusive sets of purely real and imaginary ones 
(if all kinds of dissipation are neglected). In addition, 
the crystal symmetry forces the tensor elements to de- 
couple into mutually exclusive sets of elements odd and 
even in magnetization-reveTsal, these two divisions are 
equivalent in the absence of conventional dissipation, i.e. 
real (imaginary) elements are even (odd) in the mag- 
netization. These are the results of a purely quantum- 
mechanical approach, where the Hamiltonian is Hermi- 
tian (non-dissipative) . However, the nonlinear suscepti- 
bility tensor describes the observed process of SHG, and 
thus one should not apply uniquely microscopic conclu- 
sions to the analysis of these tensor elements. Conse- 
quently, taking into account the dissipation in frequency 
space (i.e. redistribution of the response frequencies), 
will prevent the classification of tensor elements as purely 
real or imaginary ones, although for systems with higher 
symmetry the classification of tensor elements as odd and 
even ones in the magnetization (or in the antiferromag- 
netic order parameter L) can still apply jBOt . 

In summary, we have shown that the time-reversal op- 
eration, often used for the symmetry classification of 
magneto-optical phenomena, in general cannot be ap- 
plied to nonlinear optics. It should rather be replaced 
by spatial operations, resulting then in a proper descrip- 
tion of the phenomena. 
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FIG. 1. Time-reversal asymmetry in SHG. Panel (a) presents the original process, panel (b) a process in reversed time 
which would restore the symmetry, panel (c) presents a physically valid process described in reversed time. 
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FIG. 2. Light intensity distribution on the input (a) and on the output (b) of the SHG 
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FIG. 3. Nonlinear susceptibility tensor elements and resulting SHG intensity using time-reversal (panel (a)) and spin-reversal 
(panel (b)). Position of the points "A" and "B" is given by (Xe^"*)^ + ix^'')^ i 2x'^^ ■ s-nd the distance of the points "A" 
and "B" from the origin of the complex plane corresponds to the intensity of SHG from the domains A and B, respectively 
(see inset for an example of domains in Cr2 03). For simplicity, the moduli of the tensor elements have been taken as equal to 
1, but the argumentation also holds in the general case. 
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